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I.  INTRODUCTION 

Probe-corrected  near-field  techniques  have  been  widely 
used  for  the  past  40  years  to  characterize  antennas  and  trans¬ 
ducers  from  measurements  on  planar,1-5  cylindrical,6-10  and 
spherical11-20  scanning  surfaces.  The  field  of  the  antenna  or 
transducer  is  first  measured  with  a  known  probe  on  the  scan¬ 
ning  surface  in  the  near  field.  Probe-corrected  formulas  are 
then  applied  to  the  measured  near  field  to  get  the  desired  far 
field  of  the  antenna  or  transducer.  The  measurements  are 
typically  performed  in  anechoic  chambers.  The  electromag¬ 
netic  probe-corrected  formulas  have  been  derived  and  widely 
implemented  for  scanning  on  planar,  cylindrical,  and  spheri¬ 
cal  surfaces.  The  corresponding  acoustic  formulas  have  only 
been  derived  for  scanning  on  planar  and  spherical  surfaces. 

The  far  field  of  electromagnetic  and  acoustic  scatterers 
can  also  be  determined  from  near-field  measurements  in 
anechoic  chambers.  A  far-field  source  or  a  compact- 
range  reflector  can  provide  the  plane- wave  field  required  to 
illuminate  or  insonify  the  scatterers.  To  compute  the  scat¬ 
tered  far  field,  the  incident  field  and  background  interactions 
must  be  extracted  from  the  measured  total  field.  The  simplest 
method  for  eliminating  these  field  contributions  is  standard 
background  subtraction  which  involves  two  measurements. 
The  first  measurement  is  carried  out  with  the  scatterer 
present,  and  the  second  measurement  is  carried  out  with  the 
scatterer  removed.  The  results  of  the  two  measurements  are 
then  subtracted  to  get  an  approximation  for  the  scattered  near 
field,  which  can  be  processed  with  the  probe-corrected  for¬ 
mulas  to  get  the  scattered  far  field. 

It  is  appropriate  to  mention  the  related  area  of  research 
known  as  near-field  acoustical  holography,  where  near-field 
measurements  are  employed  to  backpropagate  the  pressure 
field  in  space  and  time  toward  the  source.24-27  One  can 
achieve  super-resolution  with  near-field  holography  when 
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part  of  the  evanescent  spectrum  is  captured  during  the  scan. 
In  this  paper,  however,  we  are  concerned  only  with  comput¬ 
ing  the  field  of  the  source  outside  the  scanning  surface  and 
do  not  attempt  to  backpropagate  the  near-field  data. 

In  this  paper  we  derive  the  probe-corrected  formulas  for 
cylindrical  near-field  scanning  that  can  be  used  to  compute 
acoustic  far  fields  of  transducers  and  scatterers.  We  present 
computation  schemes  and  sampling  theorems  that  allow  the 
far  fields  to  be  computed  from  sampled  values  of  the  probe 
output  on  the  scan  cylinder.  In  practice,  the  measurements 
are  carried  with  a  mechanical  scanner  over  a  truncated  scan 
cylinder  that  covers  only  a  limited  region  of  the  infinite  scan 
cylinder  that  is  used  for  the  derivation.  Therefore  we  care¬ 
fully  examine  the  probe-corrected  formulas  to  determine 
how  truncation  edges  may  affect  the  far-field  accuracy.  We 
develop  edge-correction  techniques  that  reduce  the  far-field 
errors  caused  by  truncation  and  show  that  an  angular  taper 
function  is  automatically  applied  to  the  probe  output.  The 
results  of  the  paper  are  validated  through  numerical  ex¬ 
amples  involving  a  point  source  and  a  baffled  piston  trans¬ 
ducer  probe  with  a  circular  cross  section. 

The  paper  is  organized  as  follows.  In  Sec.  II  we  derive 
the  formulas  for  cylindrical  near-field  scanning  of  acoustic 
fields  with  an  arbitrary  known  probe.  These  formulas  express 
the  field  of  the  source  region  in  terms  of  the  probe  receiving 
coefficients  and  the  probe  output  over  an  infinite  scan  cylin¬ 
der.  For  use  in  the  numerical  simulations,  we  derive  exact 
expressions  for  the  probe  output  when  the  probe  is  a  baffled 
circular  receiving  piston  transducer. 

We  consider  a  truncated  scan  cylinder  in  Sec.  El  and 
develop  edge-corrected  formulas  that  compensate  for  the 
end-point  contributions  to  the  far  field  from  the  top  and  bot¬ 
tom  edges.  We  further  study  the  angular  taper  function  that  is 
inherent  in  the  cylindrical  expansions.  In  Sec.  IV  the  far-field 
formulas  are  converted  into  discrete  computation  schemes. 
In  Sec.  V  we  consider  a  point  source  at  the  origin  to  illustrate 
the  implications  of  truncating  the  scan  cylinder.  Section  VI 
presents  conclusions. 
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FIG.  1 .  Probe  measures  the  field  of  the  source  on  the  truncated  scan  cylinder 
given  by  p=a,  - z0<z<z0 ,  <f>0<  (f>0. 


II.  THEORY  OF  CYLINDRICAL  NEAR-FIELD 
SCANNING 

We  now  derive  the  formulas  for  computing  the  acoustic 
field  of  a  source  of  finite  extent  from  either  the  pressure  field 
(Sec.  II  A)  or  from  the  output  of  a  known  probe  (Sec.  II B) 
on  a  scan  cylinder.  The  source  can  either  be  a  transducer  or  a 
scatterer.  For  scatterers  we  assume  that  background  subtrac¬ 
tion  has  been  performed,  so  that  the  incident  field  and  back¬ 
ground  reflections  have  been  removed.  The  part  of  space  not 
occupied  by  the  sources  is  a  linear  lossless  fluid  in  which  the 
acoustic  field  can  be  fully  described  by  its  pressure  field 

p(  r). 

As  shown  in  Fig.  1,  the  field  is  measured  on  the  cylinder 
p=a ,  where  p  is  the  cylindrical  coordinate  that  equals  the 
distance  to  the  z  axis.  The  scan  cylinder  in  Fig.  1  is  truncated 
and  covers  the  region  (f>^27r- cf)0.  For  the 

derivations  of  this  section  the  scan  cylinder  is  assumed  infi¬ 
nite  (zo=+°°  and  <f>\ o=0). 

We  define  the  minimum  cylinder  p=Rmin  such  that  the 
maximum  (supremum)  value  of  the  coordinate  p  for  all 
points  on  the  source  equals  /?min.  Note  that  Rmin  depends  on 
the  size  of  the  source  as  well  as  on  its  location.  For  example, 
the  value  of  Rm in  is  equal  to  V*2+y2  for  a  point  source  at 
(jc,y,z)  =  (xj ,yi,Zi).  We  assume  that  the  scan  cylinder  en¬ 
closes  the  minimum  cylinder. 

A.  The  far  field  in  terms  of  the  field  on  the  scan 
cylinder 

The  formulas  that  express  the  field  from  a  confined 
source  in  terms  of  the  pressure  on  an  enclosing  cylinder  will 
be  derived  in  this  section  from  the  standard  cylindrical  wave 


expansion.  With  e  l(at  time  dependence  suppressed,  the  basic 
outgoing  solution  p^(r)  to  the  scalar  wave  equation 


(V2  +  k2)p{r)  =  0 

(1) 

p(3Kp,<t>,z)  =  H<J\kpp)eln‘t>eikS, 

(2) 

where  k  is  the  wave  number,  kz  is  a  real  parameter,  n  is  an 
integer,  and 

kp=  Vfc2-^  (3) 

is  a  complex  parameter  with  non-negative  real  and  imaginary 
parts.  Moreover,  H^\kpp)  is  the  Hankel  function  of  the  first 
kind  and  order  n ,  and  (p,<£,z)  are  the  standard  cylindrical 
coordinates  given  in  terms  of  rectangular  coordinates  {x,y  ,z) 
as 


x  =  pcos<f),  y  =  psin<f),  z  =  z. 


(4) 


Outside  the  minimum  cylinder  of  radius  R^,  the  field  /?( r) 
of  any  source  of  finite  extent  can  be  expressed  as  a  super¬ 
position  of  the  basis  fields  in  (2)  as 


28 


+0°  /*+oo 

pip,  =  2  ein*  Fn(kz)H«\kpp)eik*dkz, 

n=— oo  J  —oo 


P  ^  ^min 


(5) 


where  Fn(kz)  is  the  spectrum  that  characterizes  the  source. 

The  formula  (5)  determines  the  field  everywhere  outside 
a  cylinder  of  radius  in  terms  of  an  integration  and  a 
summation  involving  the  spectrum  Fn(kz) .  By  applying  the 
method  of  stationary  phase,29  one  can  show  that  for  observa¬ 
tion  points  far  form  the  source  the  field  can  be  expressed  in 
spherical  coordinates  as 


p(r,0,<f>) 


2el 


ikr 


2  Fn(k  cos  e)ein4,e-Mn+l)a, 


(6) 


where  (r,  6,  <fi)  are  the  standard  spherical  coordinates  related 
v  rectangular  co  idmates  (x,y,z)  through 

x-r  sin  0cos  <f>,  y  =  r  sin  6 sin  0,  z-r  cos  0.  (7) 

Unlike  the  general  formula  (5),  the  far-field  formula  (6)  does 
not  involve  an  integration  over  kz. 

An  expression  for  the  spectrum  Fn(kz )  in  terms  of  the 
field  p(a,<f>,z)  on  the  scan  cylinder  is  obtained  by  multiply¬ 
ing  Eq.  (5)  by  e~in and  integrating  over  the  scan  cyl¬ 
inder 

^  r+°°  c2tt 


&  ^min> 

where  we  have  used  the  orthogonality  relations 


/: 


ei(kz  kz)z^z  _  27 rS(kz  -kz) 


(8) 


(9) 


and 
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FIG.  2.  Probe  situated  in  a  probe  coordinate  system.  The  input  and  output  of 
the  probe  are  determined  by  the  mode  amplitudes  ap  and  bp,  respectively. 
The  total  field  outside  the  enclosing  cylinder  for  the  probe  is  described  in 
terms  of  outgoing  and  standing  waves. 


n  =  nr 
n  =£  n' . 


(10) 


Even  though  the  radius  a  of  the  scan  cylinder  appears  on  the 
right-hand  side  of  the  expression  (8)  for  spectrum,  the  ex¬ 
pression  (8)  is  independent  of  a.  The  formulas  (5),  (6),  and 
(8)  express  the  field  everywhere  outside  a  minimum  cylinder 
that  encloses  the  sources  in  terms  of  the  field  on  the  scan 
cylinder. 


B.  Far  field  in  terms  of  the  probe  output  on  the  scan 
cylinder 

We  shall  next  present  a  simple  straightforward  deriva¬ 
tion  of  the  probe-corrected  formulas  that  express  the  field 
outside  the  scan  cylinder  in  terms  of  the  output  of  a  known 
probe  on  the  scan  cylinder.  We  shall  employ  the  outgoing 
and  standing  wave  basis  functions  of  the  source  scattering 
matrix  formulation.  The  schematic  for  the  cylindrical  near¬ 
field  measurement  system  is  shown  in  Fig.  1.  First  we  char¬ 
acterize  the  probe  with  respect  to  cylindrical  waves. 


1.  Characterization  of  the  probe 

Start  by  defining  a  probe  coordinate  system  in  which  the 
probe  is  fixed  as  shown  in  Fig.  2.  The  origin  of  this  coordi¬ 
nate  system  is  the  reference  point  of  the  probe.  The  orienta¬ 
tion  of  the  probe  with  respect  to  a  global  coordinate  system 
is  defined  in  terms  of  the  orientation  of  the  probe  coordinate 
system  with  respect  to  the  global  coordinate  system.  For  cy¬ 
lindrical  scanning  we  assume  that  the  z  axis  of  the  probe 
coordinate  system  is  parallel  to  the  z  axis  of  the  global  coor¬ 
dinate  system  that  is  used  to  describe  the  scanning  geometry. 

The  probe  is  attached  to  a  waveguide  that  supports  only 
one  propagating  mode.  The  input  to  the  probe  is  determined 
by  the  mode  amplitude  ap  and  the  output  of  the  probe  is 
determined  by  the  mode  amplitude  bp.  We  refer  to  bp  as  the 
probe  output  and  assume  that  the  probe  waveguide  is  per¬ 
fectly  matched  (ap= 0)  when  the  probe  measures  the  field  of 
the  source. 

The  field  outside  the  cylinder  in  Fig.  2  that  encloses  the 
probe  can  be  described  in  terms  of  the  outgoing  waves 


P°XP,  (f>,z)  =  (1 1) 

and  the  standing  waves 

pm{pM  =  Jn{kpp)ein^,  (12) 

where  Jn(kpp)  is  the  Bessel  function  of  order  n.  The  probe 
receiving  coefficients  can  now  be  defined  as  follows.  If  the 
total  field  outside  the  enclosing  cylinder  is  the  single  stand¬ 
ing  cylindrical  wave  Jn(kpp)ein<f>elkzZ ,  the  probe  output  is 
given  by  bp=Cn(kz ),  where  Cn(kz)  are  the  probe  receiving 
coefficients. 

The  probe  receiving  coefficients  will  now  be  expressed 
in  terms  of  the  probe’s  plane-wave  receiving  characteristic 
Rp(kx,ky ),  defined  to  be  the  probe  output  when  the  incident 
field  is  the  plane-wave  el^kxX+kyy+kzZ\  with  k2x+k2+k2-k2.  The 
plane-wave  receiving  characteristic,  which  is  also  defined 
with  the  probe  located  in  the  probe  coordinate  system,  can  be 
expressed  both  in  terms  of  the  plane-wave  transmitting  spec¬ 
trum  and  the  far-field  pattern  of  the  probe  when  the  probe  is 
a  reciprocal  electroacoustic  transducer.30,5  For  example,  if 
the  probe  is  reciprocal  and  its  far  field  is  expressed  in  terms 
of  a  far-field  pattern  Tp{Q,  <£>)  as  p(r)  ~  Tp{ 0,  cf))elkrf  r,  the 
receiving  characteristic  is 


Rp(k  cos  (f>  sin  0,k  sin  (f>  sin  9)  =  - 


Tp(lt-  (9,77+  <f>) 

ivYpPmo 


(13) 


where  Yp  is  the  characteristic  admittance  for  the  propagating 
mode  of  the  probe  waveguide  feed  and  pm0  is  the  mass  den¬ 
sity  of  the  undisturbed  fluid. 

With  the  angles  <j)k  and  0k  defined  such  that  kx 
=kpcos  <f>h  ky=kp sin  <f>h  and  kz=k  cos  0h  the  cylindrical  ex¬ 
pansion  of  the  incident  plane-wave  el^kxX+kyy+k^  is 
+00 

ewv+kyy+k7z)=  2  I nJn(kpP)ein^Wc*.  (14) 

n=- oo 


Note  that  kx< 0  and  7r/2<(f>k<3v/2  for  plane-wave  com¬ 
ponents  generated  by  sources  that  are  located  in  the  half 
space  x>0.  By  inspection,  the  plane- wave  receiving  charac¬ 
teristic  is  found  to  be 

+oo 

Rp(kx,ky)=  2  inCn(kz)e-in^,  (15) 


which  in  turn  yields  the  desired  expression  for  the  receiving 
coefficients 


rn  r 2 * 

:"(kz)  =  ^J0  1 


Rp{k  cos  <f>k,k  sin  <f)k)ein^d<f>k.  (16) 


As  will  be  demonstrated  below,  the  far  field  of  the  source 
can  be  determined  from  the  plane- wave  receiving  character¬ 
istic  in  the  region  of  propagating  waves  where  kx,  ky,  and  kz 
are  real. 


2 .  Piston  probe 

We  consider  now  a  baffled  circular  receiving  piston 
transducer  whose  piston  resides  in  the  y-z  plane  of  the  probe 
coordinate  system.  For  brevity  we  shall  refer  to  this  receiv- 
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FIG.  5.  Geometry  for  translation  and 
rotation  of  cylindrical  wave  functions. 


3.  Probe-corrected  formulas 

The  probe-corrected  formulas  for  the  field  of  the  source 
will  be  derived  in  this  section  using  a  straightforward  ap¬ 
proach  that  involves  two  steps.  First,  an  expression  is  derived 
for  the  probe  output  bp(a,<f>,z )  in  terms  of  the  unknown 
expansion  coefficients  Fm(kz )  for  the  source  and  the  known 
probe  receiving  coefficients  Cn(k z).  Second,  this  expression 
for  the  probe  output  is  inverted  to  get  the  expansion  coeffi¬ 
cients  of  the  source  in  terms  of  the  probe  receiving  coeffi¬ 
cients  and  the  probe  output. 

Multiple  interactions  between  the  source  and  the  probe 
will  be  neglected.  That  is,  the  field  scattered  by  the  probe, 
rescattered  by  the  source,  and  returned  to  the  probe,  does  not 
change  the  output  of  the  probe.  The  effect  of  multiple  inter¬ 
actions  can  be  taken  into  account  formally  by  the  source 
scattering  matrix  formulation.8  However,  this  formulation 
does  not  provide  quantitative  information  about  these  mul¬ 
tiple  interactions,  and  to  obtain  useful  probe-corrected  for¬ 
mulas  they  have  to  be  neglected. 

To  derive  an  expression  for  the  probe  output  when  the 
probe  is  illuminated  by  the  field  of  the  source,  we  need  to 
express  the  cylindrical  waves  in  Eq.  (5)  in  terms  of  cylindri¬ 
cal  waves  defined  in  the  probe  coordinate  system.  To  distin¬ 
guish  between  the  global  coordinate  system  (in  which  the 
scan  cylinder  is  given  by  p=a)  and  the  probe  coordinate 
system  (in  which  the  probe  is  fixed)  we  attach  primes  to  all 
coordinates  defined  with  respect  to  the  probe  coordinate  sys¬ 
tem. 

Assume  that  the  reference  point  of  the  probe  is  at 
( a,<f)s,zs )  on  the  scan  cylinder.  Then  we  need  to  transform 
the  cylindrical  wave  functions  in  Eq.  (5)  into  cylindrical 
wave  functions  in  the  probe  coordinate  system,  which  is  cen¬ 
tered  at  ( a,(f>s,zs )  with  its  x  axis  pointing  toward  the  center 
of  the  global  coordinate  system.  To  achieve  this  transforma¬ 
tion  we  use  the  rotation-translation  formula  for  cylindrical 
wave  functions 


+00 

=  (-  ly V**  2  HWm(ps)Jn(p')ein+',  ft>p\ 

n=-o o 

(25) 

where  the  primed  quantities  are  defined  in  the  probe  coordi¬ 
nate  system  and  the  unprimed  quantities  are  defined  in  the 
global  coordinate  system,  as  shown  in  Fig.  5. 

Equation  (25)  can  be  inserted  into  expression  (5)  to  get 
the  following  formula  for  the  pressure  on  the  scan  cylinder: 

+00 

p(a,<f>s,zs)  =  2  (-DV^* 

m=-°° 

+0°  r+oo 

X  2  [Jn(kpP>in*'eik*']Fm(kz)H«}m 
n—— 00  *'  —oo 

X(kpa)  eLK*sdkz ,  (26) 

which  expresses  the  outgoing  waves  that  emanate  from  the 
source  in  the  global  coordinate  system  in  terms  of  standing 
waves  Jn(kppr)ein<f>  einkzZ  in  the  probe  coordinate  system. 
From  the  definition  of  the  probe  receiving  coefficients,  we 
see  that  the  total  probe  output  is 

+00  +00 

bp(a,<f>,z)  =  2  (-1)V"^2  Cn(kz)Fm(kz)H(nl}m 

m~— 00  n~— 00  —oo 

X{kffl)eik^sdkz.  (27) 

The  expression  (27)  is  easily  inverted  by  use  of  the  orthogo¬ 
nality  relations  (9)  and  (10)  to  get  the  final  probe-corrected 
expression  for  the  spectrum 

Fm(kz)  =  Dm(kz)Im{kz),  (28) 

where 
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FIG.  6.  Correction  factor  at  kz=0  for  piston  probes  with  h= 0.1X,  /i=0.61X,  h=k,  and  h=2\  that  measure  the  field  on  a  scan  cylinder  of  radius  a- 30X. 


r+ oo  r2ir 

Im(kz)  =  I  bp(a,(f>s,zs)e~,m^e~,k^d(t>sdzs  (29) 
J-  oo  JO 

is  a  Fourier  integral  of  the  probe  output  over  the  scan  cylin¬ 
der  and 

(-  \)m 

Dm{kz)  = - ^  -  (30) 

4^2  Cn(kz)H^lm(kpa) 

n=-oo 

is  a  correction  factor  that  depends  on  the  probe  receiving 
coefficients  and  the  scan  cylinder  radius.  The  far  field  of  the 
source  is  expressed  in  Eq.  (6)  in  terms  of  the  expansion 
coefficients  Fm(kz)  evaluated  at  kz-k  cos  0 ,  and  the  formula 
(16)  determines  the  probe  receiving  coefficients  Cn(kz)  in 
terms  of  the  plane-wave  receiving  characteristic  of  the 
probe.  For  an  ideal  probe  whose  output  equals  the  pres¬ 
sure  (Cn=0  for  ni=0  and  C0=  1),  Eq.  (28)  reduces  cor¬ 
rectly  to  the  nonprobe-corrected  formula  (8). 

Figure  6  shows  the  magnitude  in  decebels  of  the  correc¬ 
tion  factor  (30)  for  piston  probes  with  ft=0.1\,  /*=0.61\,  h 
=\,  and  h- 2\.  The  scan  cylinder  radius  is  a- 30A  and  the 
observation  direction  is  6-  90°,  corresponding  to  kz=0.  We 
shall  now  describe  how  the  correction  factor  corrects  for  the 
probe  pattern  for  far-held  observation  points  near  0=7t/2. 

First  assume  that  the  field  is  produced  by  a  point  source 
at  the  origin  so  that  /?min=0.  Then  the  expansion  coefficients 
Fm( 0)  are  zero  for  0,  and  only  the  correction  factor 
Z>0(0)  comes  into  play.  For  points  on  the  scan  cylinder  near 
the  x-y  plane,  the  output  is  roughly  the  same  for  all  four 
probes  because  the  point  source  is  in  the  direction  of  the 
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main  beam  (see  Fig.  3  with  /3=0).  The  far  field  of  the  point 
source  at  6=  tt12  is  mainly  determined  by  the  probe  output 
near  the  x-y  plane,  so  the  correction  factor  D0( 0)  should  be 
roughly  the  same  for  the  four  probes,  as  confirmed  by  Fig.  6. 

Next  assume  that  the  field  is  produced  by  a  point  source 
on  the  x  axis  at  x=3a/4,  for  which  the  expansion  coefficients 
Fm( 0)  with  \m\  <  150  are  non- negligible.  The  probe  sees  this 
point  source  at  an  angle  /3=37°  when  the  probe  is  at 
(x,y,z)=(0,<z,0)  and  (x,y,z)  =  (0,-«,0).  At  all  other  obser¬ 
vation  points  in  the  x-y  plane,  the  probes  sees  the  point  scat- 
terer  at  an  angle  (3  that  is  less  than  37°.  Consider  the  piston 
probe  with  h=k.  From  Fig.  (3)  we  see  that  (3=31°  is  near  the 
first  null  of  this  probe,  so  very  low  outputs  are  obtained 
when  the  probe  is  near  (x,y ,z)  =  (0,a,0)  and  (x,y,z)  =  (0, 
-a, 0).  To  compensate  for  these  low  outputs,  the  correction 
factor  Dm( 0)  has  a  peak  at  \m\- 120  and  attains  large  values 
around  this  peak.  Hence,  to  compute  some  of  the  higher- 
order  expansion  coefficients  of  the  source  Fm( 0),  the  integral 
over  the  scan  cylinder  Im{ 0)  must  be  boosted  by  a  large  cor¬ 
rection  factor  Dm( 0).  These  arguments  can  be  repeated  for 
the  other  piston  probes  in  Fig.  3. 

This  discussion  shows  that  nulls  and  low  sidelobe  levels 
in  the  probe  receiving  pattern  result  in  large  values  of  the 
correction  factor.  Correction  factors  that  vary  significantly 
with  m,  can  often  lead  to  errors  in  the  computed  far  field 
because  some  of  the  higher-order  expansion  coefficients  of 
the  source  Fm  are  computed  as  the  product  of  a  small  number 
that  can  be  relatively  inaccurate  (the  integral  of  the  probe 
output  over  the  scan  cylinder  Im)  and  a  very  large  number 
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(the  correction  factor  Dm).  Therefore,  it  is  advantageous  to 
avoid  probes  that  have  nulls  and  low  sidelobe  levels  in  the 
region  of  the  sources. 

In  Sec.  IV  we  show  how  the  far  field  can  be  computed 
from  sampled  values  of  the  probe  output  on  the  scan  cylin¬ 
der.  If  the  angular  sampling  rate  is  A</>,  the  far-field  summa¬ 
tion  (6)  includes  the  terms  with  \n  \  <  ir/Acf?  (recall  that  the 
spectrum  in  Eq.  (6)  is  given  in  Eq.  (28)  as  Fn=DnIn).  Hence, 
by  decreasing  the  sample  spacing  A  </>  the  far-field  error  may 


increase  because  more  of  the  peaks  of  the  correction  factor 
are  included.  This  phenomenon  is  illustrated  in  Sec.  V.  Of 
course,  the  sample  rate  A </>  must  always  be  chosen  small 
enough  to  properly  resolve  the  field  of  the  source. 

We  shall  now  briefly  discuss  an  approximate  probe- 
corrected  formula  obtained  by  inserting35 

rfn-Jkpd)  =*  (-  l)m(-  i)ne~in  arcsin(my[V»//m1)(V)  (3  0 
into  Eq.  (30)  to  get 


Dm(kz)  - 


_ 1 _ 

4'7T2Rp(k,pcos[arcsin(m/[kpa])  +  7r],/:psin[arcsin(m/[£pa])  +  ir])fl<m1)(fepa) 


(32) 


An  electromagnetic  analog  of  Eq.  (32)  was  derived  by 
Borgiotti9,10  and  subsequently  rederived  by  Yaghjian.8  The 
approximate  formula  (32)  expresses  the  correction  factor 
Dm(kz)  directly  in  terms  of  the  receiving  characteristic 
Rp(kx,ky)  and  avoids  the  use  of  the  cylindrical  receiving  co¬ 
efficients  Cn(kz )  that  occur  in  the  exact  formula  (30).  By 
invoking  the  relation  (13),  the  formula  (32)  expresses  the 
correction  factor  directly  in  terms  of  the  probe  far-field  pat¬ 
tern.  For  kpa>\,  expression  (31)  is  accurate  only  for  small 
\n\  and  values  of  \m\  that  are  not  too  close  to  kpa.  Hence,  Eq. 
(32)  is  accurate  only  for  small  probes  with  only  a  few  non¬ 
zero  receiving  coefficients  Cn(kz)  and  for  sources  with  R^n 
not  too  close  to  a.  The  approximate  formula  (32)  is  sin¬ 
gular  when  arcsin(m/[&pfl])  +  7r  coincides  with  a  null  of  the 
receiving  characteristic.  These  singularities  correspond  to 
the  “peaks”  discussed  above  where  the  exact  correction 
factor  attains  large  values  but  remains  finite.  The  approxi¬ 
mate  formula  (32)  thus  fails  to  correctly  predict  the  values 
of  the  correction  factor  near  its  peaks.  However,  it  is  a 
valuable  formula  that  avoids  the  use  of  the  cylindrical 
receiving  coefficients  and  produces  accurate  far  fields  for 
small  probes  whose  patterns  are  null  free.  As  we  shall 
demonstrate  below,  small  probes  have  additional  advan¬ 
tages. 

The  fact  that  the  exact  correction  factor  is  without  sin¬ 
gularities  distinguishes  cylindrical  scanning  from  planar 
scanning:  In  planar  scanning  the  exact  correction  factor  is 
proportional  to  1  / Rp  and  thus  has  singularities  at  the  nulls  of 
Rp,  as  proven  by  Hansen  and  Yaghjian5.  This  discussion  was 
confined  to  far-field  observation  directions  near  0=7t/2.  In 
Sec.  V  we  shall  show  how  nulls  in  the  probe  pattern  affects 
the  computed  far  field  for  observation  points  away  from  6 
=  tt/2. 

III.  TRUNCATION  OF  THE  SCAN  CYLINDER 

The  exact  expressions  (8)  and  (28)-(30)  for  the  expan¬ 
sion  coefficients  Fm(kz)  involve  an  integration  over  an  entire 
infinite  scan  cylinder  of  either  the  pressure  or  the  probe  out¬ 
put.  Practical  near-field  measurements  are  carried  out  only 
over  a  truncated  scan  cylinder,  which  we  assume  is  given  by 


-Zo ^Z^Zo  and  where  z0>0  and  0^ </>0 

<  7 r.  The  near-field  data  required  by  Eqs.  (8)  and  (28)-(30)  is 
therefore  not  available  in  practice.  Fortunately,  experience 
has  shown  that  very  accurate  far  fields  can  be  obtained  in 
certain  regions  from  measurements  on  a  truncated  scan  cyl¬ 
inder.  In  this  section  we  shall  study  the  errors  caused  by 
truncation  and  develop  methods  for  reducing  them.  Trunca¬ 
tion  in  the  z  and  <j>  directions  are  considered  in  Secs.  HI  A 
and  HI  B,  respectively. 

First  we  state  a  result  from  Hansen  et  al.22  that  restricts 
the  range  of  observation  angles  for  which  we  can  accurately 
compute  the  far  field  from  a  truncated  scan:  If  the  scanned 
area  is  thought  of  as  a  transparent  surface  and  the  rest  of  the 
infinite  scan  cylinder  is  opaque ,  then  accurate  far  fields  can 
be  achieved  only  in  directions  from  which  one  can  see  the 
entire  source  region.  Similar  results  have  been  cited  for  pla¬ 
nar  and  spherical  scanning  surfaces.36-39 


A.  Truncation  in  the  z  direction 

Typically  the  field  outside  the  scanned  area  is  neglected 
and  the  exact  expression  (29)  is  approximated  by 

fz o  r27T-0o 

/.(*«)“  bp{a,<t>s,zs)e-in+°e-ik*d<l>4zs.  { 33) 

J  -Zo  J  </>o 

It  is  thus  assumed  that  the  contribution  to  the  integral  (29) 
from  the  region  not  covered  by  the  scan  is  negligible.  For 
cylindrical  near-field  measurements  of  directive  transducers 
this  assumption  is  often  valid  because  the  measurements  are 
set  up  so  that  very  little  energy  is  radiated  in  directions  not 
covered  by  the  truncated  scan  cylinder.  However,  a  nondi¬ 
rective  transducer  or  a  scatterer  may  radiate  significantly  in 
all  directions,  including  the  ones  that  are  not  covered  by  the 
scan  cylinder.  Hence,  it  is  not  always  a  valid  approximation 
to  neglect  the  contribution  to  the  integral  in  Eq.  (29)  from  the 
region  not  covered  by  the  scan. 

We  shall  next  present  two  sets  of  asymptotic  correction 
formulas  that  approximate  the  contribution  from  the  region 
of  the  scan  cylinder  with  \z  \  >Zq.  Detailed  derivations  can  be 
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found  in  Hansen  et  al.  First  we  assume  that  the  field  on  the 
scan  cylinder  behaves  roughly  as  a  plane  wave  elk ^  as  z 
-+  ±°°  and  find  that 


n2TT-<f>0 

bp{a,  <t>s,zs)e~m^e~lk^d<l)sdzs 

br\ 


fZo  r2TT-<f>0 

-z0  J  <f> 0 

eikzz  0  r'Z'K-i 

i{k  +  &z)  J_0o 

e-ikzzo  rhr-4 

Kk-kJJ, 


bp(a,</>s,- z0)e  m+*d<f>. 


bp(a,<f>s,zo)e  m^d^s.  (34) 


We  use  the  term  “edge-corrected”  formula  because  the  inte¬ 
gral  over  the  uncovered  areas  of  the  infinite  scan  cylinder  are 
approximated  by  integrations  along  the  top  and  bottom  edges 
of  the  truncated  scan  cylinder. 

The  second  set  of  correction  formulas  assumes  that  field 
behaves  as  a  spherical  wave  elkr  as  z— >  ±°°.  The  edge- 
corrected  formula  for  the  spherical-wave  assumption  is40 


■z0  r2n-4>0 

bp(a,  <f>s,zs)e~in<l’se~ik^sd(f>sdzs 

~Z(\  J  <bn 

eikzZ0 

r2rr-(f)Q 

|  bp(a,<f>s,-zo)e~in4’sd<f>s 

>0 

i(k  cos  0O  +  k:) . 

e~ikzz0  | 

r2ir-<f>o 

1  bp(a,  <f>s,zo)e~m<t’sd<f>s. 

i(k  cos  0O  -  kz )  J 

(35) 


where  cos  0o=zo/  ^1  a2 +z%.  Note  that  the  edge-corrected  for¬ 
mula  (35)  is  singular  for  far-field  observation  directions 
with  0=  0q  and  0=7r/2-0o.  This  follows  from  the  far-field 
expression  (6)  which  shows  that  the  far  field  in  the  direc¬ 
tion  ( 0, fi)  is  determined  from  the  spectrum  evaluated  at 
kz=k  cos  0.  The  edge-corrected  formula  (34)  does  not  have 
these  singularities. 


B.  Truncation  in  the  fi  direction 

For  the  angular  truncation  of  the  scan  cylinder  we  can¬ 
not  derive  asymptotic  correction  formulas  similar  to  Eqs. 
(34)  and  (35)  because  no  asymptotic  expression  suggests  it¬ 
self  for  the  fi  dependence  of  the  field.  Indeed,  the  fi  depen¬ 
dence  of  the  field  on  the  scan  cylinder  depends  strongly  on 
the  source  shape,  so  we  do  not  have  a  general  expression  for 
bp(a,fi,z)  that  holds  when  fi  is  outside  the  scanned  region 
0o^0^27r“0o- 

To  determine  the  effect  of  0  truncation,  insert  the  for¬ 
mula  (28)  for  the  spectrum  Fm(kz )  into  the  far-field  formula 
(6)  to  get 


p(r,0,0) 


eikr  r+°°  r2ir 

r  oo  Jn 


bp{a,  cf>s,zs)G(k,a,  6,  <f> 


-  &;{£„})*' 


,-ikzs cos  0 


d(f>sdzs  > 


(36) 


where 


+0°  /_  g-i{m+l)Tr/2 ^inup 

G(k,a,0,<f>\{Cn})  =  2  - ^—L - 

2  u2  2  Cn(k  cos  9)H^}m(ka  sin  0) 

n=- oo 

(37) 

will  be  called  the  taper  function.  Equation  (36)  expresses  the 
far  field  of  the  source  directly  in  terms  of  an  integration  over 
the  scan  cylinder  of  the  probe  output  multiplied  by  the  taper 
function.  This  expression  is  useful  for  determining  the 
0-truncation  errors  but  not  efficient  for  actually  computing 
the  far  field  since  the  fast  Fourier  transform  (FFT)  cannot  be 
fully  exploited.  We  shall  now  examine  the  taper  function 
(37)  for  an  ideal  probe  and  for  piston  probes. 


1.  Ideal  probe 


For  an  ideal  probe  (Cn= 0  for  n  ¥=  0  and  C0=  1 ),  the  taper 
function  (37)  reduces  to 


Gfka ,  0, 0) 


+00 

2 


m=- oo 


g-i(m+ 1 )  Trllgim  <f> 

iT^H^Xka  sin  0) 


(38) 


The  piston  probe  with  /i=0.1\  is  approximately  an  ideal 
probe.  By  inspection,  we  see  that  the  function 
Gi(ka,0,(j)-(l)^e~lkZsCOS  0  equals  the  total  field  on  a  hard 
circular  cylinder  that  is  illuminated  by  a  plane  wave  with 
amplitude  (27 ri)~lka  sin  0  that  propagates  in  the  -f  direc¬ 
tion.  Here,  r  is  the  unit  vector  for  the  far-field  observation 
direction.  Hence,  the  far  field  of  the  source  is  obtained  by 
integrating  over  the  scan  cylinder  the  product  of  the  near 
field  of  the  source  and  the  total  field  on  a  hard  cylinder , 
which  coincides  with  the  scan  cylinder  and  is  illuminated  by 
a  plane  wave  that  originates  at  the  far-field  observation 
point. 

Jones41  used  the  Watson  transform  to  prove  that  as 
ka  sin  0— *  oo  the  function  Gfka^O,  0)  decays  exponentially 
for  7t/2<|0|^77,  has  a  transition  zone  around  |0|=7t/2, 
and  behaves  as 


Gfka ,  0, 0) ' 


ka  sin  0 
2  m 


cos  fie 


,-ika  sin  0  cos  <f) 


\fi\  <  77/2,  ka  sin  0— >  +  oo  .  (39) 

Figure  7  shows  Gt{ka  ,0,0)  and  its  asymptotic  approxima¬ 
tion  (39)  for  ka  sin  0=25.  The  asymptotic  approximation 
works  well  even  at  this  low  value  of  ka  sin  0. 

The  qualitative  behavior  of  Gt(ka,  0,  fi)  plotted  in  Fig.  7 
for  ka  sin  6=25  remains  the  same  for  all  larger  ka  sin  0. 
Hence,  the  far-field  formula  (36)  shows  that  the  field  in  a 
given  direction  with  angle  fi  is  determined  by  integrating  the 
field  on  the  scan  cylinder  multiplied  by  an  oscillating  func¬ 
tion  of  fis  that  is  symmetric  and  decays  away  from  fis=  fi.  In 
other  words,  the  cylindrical  wave-function  expansion  of  the 
far  field  automatically  applies  an  angular  taper  function  to 
the  field  on  the  scan  cylinder.  No  taper  function  is  applied  for 
the  zs  dependence.  The  asymptotic  expression  (39)  shows 
that  the  oscillations  increase  in  frequency  with  increasing 
ka  sin  0. 
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FIG.  7.  Real  (top)  and  imaginary  (bottom)  parts  of  the  function  Gt  and  its  asymptotic  approximation  for  ka  sin  0=25. 


2.  Piston  probe 

Figure  8  shows  an  upper  envelope  of  the  taper  function 
in  Eq.  (37)  for  piston  probes  with  /i=0.1\,  h- 0.61X,  h-\ , 
and  h=2\.  The  scan  cylinder  radius  is  a=30X  and  the  obser¬ 
vation  direction  is  #=90°,  corresponding  to  kz=0.  To  visual¬ 


ize  the  rapidly  oscillating  taper  function,  we  computed  an 
upper  envelope  that  equals  the  maximum  magnitude  of  the 
taper  function  over  5°  intervals. 

Let  us  first  discuss  the  behavior  of  the  taper  function  for 
/i=0.1\  (this  piston  probe  is  almost  an  ideal  probe  whose 


♦ 


FIG.  8.  Upper  envelope  of  the  taper  functions  for  a=30\  and  0=90°  for  piston  probes  with  h= 0.1\,  /z=0.61\,  h=\,  and  h-2k.  All  non-negligible  terms  of 
the  expansion  for  the  taper  function  are  included. 
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output  equals  the  pressure).  In  the  region  from  (f>=0°  to  0 
=  90°  the  taper  function  decays  at  a  moderate  rate  from  30  to 
about  5  dB.  Beyond  <f>=90°  the  decay  is  exponential 
(creeping- wave  region),  and  at  <f>=  180°  the  magnitude  of  the 
taper  function  is  below  -100  dB.  Hence,  the  probe  output  on 
the  scan  cylinder  at  points  that  lie  more  than  90°  away  from 
the  far-held  observation  direction  are  attenuated  by  more 
than  25  dB  when  the  fields  are  measured  with  a  piston  probe 
of  radius  /i=0.1X.  We  shall  later  see  that  the  strong  taper 
helps  reduce  truncation  effects. 

The  taper  functions  for  the  piston  probes  with  h 
=0.6  IX,  h=\,  and  h= 2\  do  not  exhibit  such  rapid  decay.  For 
example,  the  taper  function  for  the  piston  probe  of  radius  h 
=0.6  IX  decays  only  5  dB  over  the  entire  range  from  cf>=Q° 
to  <£=  180°,  and  thus  only  a  mild  taper  is  applied  to  the  probe 
output.  The  values  of  the  taper  functions  in  Fig.  8  were  com¬ 
puted  from  the  summation  (37)  with  all  non-negligible  terms 
included.  For  practical  far-held  calculations,  the  summation 
in  Eq.  (37)  is  truncated  at  a  mode  number  determined  by  the 
size  of  the  source  region.  This  truncation  of  the  summation 
strongly  affects  the  taper  function,  as  will  be  demonstrated  in 
Sec.  V. 

C.  Green’s  function  representation 

Let  us  briefly  discuss  a  far-held  formula  that  is  based  on 
the  free-space  Green’s  function  rather  than  the  cylindrical 
wave  functions.  From  Hansen  et  al.5  it  follows  that  the  far 
held  can  be  expressed  in  terms  of  the  pressure  and  its  normal 
derivative  on  the  scan  cylinder  as 

iaelkr  f +0°  f2ir  T 

+  <■-"''<(<// <ft’.  (40) 

dp  J 

The  normal  derivative  of  the  pressure  can  be  related  to  the 
velocity  by  the  formula  d/  dp'  p(r')=ia)pm0p'  •  v(r'),  where 
pm0  is  the  mass  density  of  the  undisturbed  fluid. 

This  formula  is  seldom  used  in  near-held  scanning  be¬ 
cause  it  requires  that  both  the  pressure  and  its  normal  deriva¬ 
tive  be  measured.  Unlike  the  formulas  that  are  based  on  the 
cylindrical  wave  functions,  the  formula  (40)  does  not  apply 
an  angular  taper  function  to  the  held  on  the  scan  cylinder. 
Hence,  angular  truncations  of  the  scan  cylinder  would  in 
general  affect  Eq.  (40)  more  than  they  affect  formulas  based 
on  cylindrical  wave  functions. 

IV.  SAMPLING  THEOREMS  AND  COMPUTATION 
SCHEMES 

The  formulas  derived  so  far  for  the  far  held  of  the  source 
involve  integrations  of  the  probe  output  over  the  scan  cylin¬ 
der  and  inhnite  summations  over  angular  modes.  In  this  sec¬ 
tion  we  discretize  these  formulas  so  that  the  far  held  of  the 
source  can  be  computed  from  a  hnite  number  of  probe- 
output  values  sampled  at  regular  intervals  on  the  scan  cylin¬ 
der.  As  part  of  the  discretization,  we  present  sampling  theo¬ 
rems  that  determine  the  sample  spacing  required  for  accurate 
far-held  calculation. 


The  scan  cylinder  is  truncated  and  given  in  cylindrical 
coordinates  by  p=a  and  |z|  ^Zo,  and  the  held  on  the  scan 
cylinder  is  measured  at  the  grid  points  expressed  in  terms  of 
the  cylindrical  coordinates  as 

<£=(€- 1)A  cf>,  €=1,2,...,A^,  (41) 

and 

z  =  (m-  l)Az-Zo>  m-  1,2, ...  ,NZ.  (42) 

The  formulas  of  this  section  require  the  probe  over  the  entire 
360°  range.  If  the  scan  cylinder  is  angularly  truncated,  one 
simply  inserts  zeros  into  the  array  positions  that  corre¬ 
spond  to  measurement  points  that  are  outside  the  scanned 
area.  Note  that  Zo  satisfies  the  relation  zo=(Az-l)Az/2. 

If  every  point  on  the  scan  cylinder  is  at  least  one  wave¬ 
length  or  so  away  from  the  nearest  point  on  a  nonresonant 
source,42  the  sampling  intervals  A  (j)  and  Az  can  be  chosen  as 
follows:17 


A^=— ,  ^  =  2int(fc/?min  +  «1), 

(43) 

AZ=? 

(44) 

where  R^n  is  the  radius  of  the  minimum  cylinder,  nx  is  a 
small  integer,  and  X  is  the  wavelength.  In  Eq.  (43)  the 
function  “int”  denotes  the  integer  value. 

The  sampling  theorems  show  that  the  infinite  summation 
(6)  can  be  replaced  by  the  finite  summation 

2  eikr  N° 

p(r ,  0,4>) - 2  Dn(k  cos  9)In(k  cos  0)em<^“*7r*n+1)/2, 


where  Dn(k  cos  9)  is  the  correction  factor  defined  in  Eq. 
(30),  and  N0  is  an  integer  that  must  be  large  enough  to 
include  all  the  significient  modes  of  a  give  source  region. 
When  strong  edge  effects  are  present,  N0  should  be  large 
enough  to  prevent  aliasing  caused  by  the  discontinuity  of 
the  near-field  data  at  the  truncation  edges.22  Since  only 
the  values  of  In(k  cos  0)  with  \n  \  ^N^/2  can  be  computed 
accurately,  one  must  chose  N0^N(j,/2. 

Next  we  introduce  the  two-dimensional  array 

bp((,m)  =  bp(a,(t-  1)A <£,(m  -  l)Az -  z0), 

€=1,2,...,AT0,  m=  1,2, ...  ,NZ,  (46) 

which  contains  the  probe  output  bp  at  the  grid  points.  A 
straightforward  discretization  of  the  formula  (33)  for  the 
spectrum  then  gives 

^  Nz 

In(kz)  =  2  e-in^(-l)e-ikMm-Vbp{t,m). 

t-\  m- 1 

(47) 

Note  that  both  the  far-field  formula  (45)  and  the  probe- 
corrected  formula  (47)  for  the  spectrum  can  be  computed 
efficiently  with  the  FFT. 
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We  can  also  discretize  the  edge-corrected  formulas  (34) 
and  (35)  to  get 

**  Nz 

In(kz)  =  AzA(f>eikzz«'Z  2 

€=1  m- 1 


2  €=1 

N, 


2  e~in^e~l^bp(£  ,NZ) 


1) 

*(fc  +  ^)€=1 

*(£-£;)  €=1  ^ 
when  the  field  is  assumed  to  behave  as  and 

/„0y  =  AzAtfe^E  2 


(48) 


AzA<f>elk^° 


=  1  m=l 


2 


f=l 


AzAtfe*W>§ 


2  e-in^(~l)bp{i,Nz) 


€=1 


A(f>eikzZo 
i(k  cos  0O  +  &z) 


2  ^a^“1)L(€,1) 


z/  e=i 


A^'^o  § 


/(A:  cos  0O  -  *z)  ^=i 


2  e-«‘nA^“i)y€,7Vz)  (49) 


when  the  field  is  assumed  to  behave  as  and  cos  0O 
=Zo/V<22+Zo-  The  second  and  third  in  Eqs.  (48)  and  (49) 
ensure  that  the  contributions  from  the  top  and  bottom 
edges  are  included  only  once.  The  number  of  operations 
required  to  compute  the  spectrum  from  any  one  of  the 
formulas  (47)-(49)  is  on  the  order  of  N<^z\og2(N<f>9Nz). 
Hence,  the  edge-correction  terms  do  not  significantly  add 
to  the  computational  effort. 


V.  FAR-FIELD  ERRORS  FOR  A  POINT  SOURCE 

In  this  section  we  shall  compute  the  far  field  of  an 
acoustic  point  source  from  near-field  measurements  on  a 
truncated  scan  cylinder.  The  field  is  measured  with  piston 
probes  of  varying  size  and  the  probe  output  is  computed 
from  the  exact  formula  (22).  Hence,  if  the  scan  cylinder  is 
untruncated,  the  computed  far  field  equals  the  exact  far  field, 
regardless  of  the  size  of  the  piston  probe.  The  acoustic  field 
is  generated  by  a  single  point  source  (i.e.,  a  nondirective 
source)  located  inside  the  scan  cylinder,  so  the  truncation 
edges  of  the  scan  cylinder  are  illuminated  by  relatively 
strong  fields,  and  the  truncation  effects  will  be  clearly  vis¬ 
ible.  The  scan  cylinder  is  given  by  a=30\  and  -40k  <z 


TABLE  I. 


Average  far-field  error  over  the  region  0 

4><2ir  for  O-ttH 

h= 0.1\  h-- 

=0.61\ 

h=k 

h=2\ 

(%) 

w 

(%) 

(%) 

No  edge  correction 

3.1 

0.7 

0.5 

0.3 

Plane-wage  edge  correction 

0.9 

0.2 

0.1 

0.1 

Spherical-wave  edge  correction 

1.3 

0.2 

0.2 

0.1 

<39.5\.  Throughout,  the  sampling  spacing  in  z  is  Az=\/2. 
First  assume  that  the  point  source  is  located  at  r5=12ky 
-5\z. 

We  begin  by  investigating  the  improvements  in  far-field 
accuracy  achieved  with  the  formulas  that  correct  for  the  trun¬ 
cations  at  the  top  and  bottom  edges  of  the  scan  cylinder.  We 
let  the  scan  cylinder  be  untruncated  in  $  and  chose 
=360  corresponding  to  an  angular  sample  spacing  of  A  <p 
=  1°.  (With  ^min=  12\,  the  required  sample  spacing  is  about 
2.3°.)  The  sample  spacing  in  z  is  Az=\/2.  We  shall  compare 
the  far  fields  obtained  with  four  different  probe  sizes:  h 
=0.1k,  /i=0.61\,  /i=\,  and  h-2X.  The  exact  far  field  of  a 
single  point  source  located  at  rs  is 


eikr 

p0(r)  ~  e-iki  T*—.  (50) 

r 

Table  I  shows  the  average  error  of  the  far  field  for  0-  it  1 2 
and  0^(f)<2ir  computed  from  this  z-truncated  scan  using 
Eqs.  (47)-(49). 

We  see  that  the  edge-correction  formulas  significantly 
reduce  the  far-field  error  for  all  sizes  of  the  piston  probe. 
Moreover,  the  error  is  smaller  for  directive  probes  because  of 
the  reduced  magnitude  of  the  probe  output  near  the  top  and 
bottom  edges  of  the  scan  cylinder.  In  the  remainder  of  this 
section  we  employ  the  edge-corrected  formula  (48)  based  on 
the  plane- wave  assumption. 

Figure  9  shows  the  far-field  error  as  a  function  of  0  for 
(f>=  0  with  the  scan  cylinder  still  untruncated  in  <fi.  The  error 
is  calculated  at  the  angles  0=arccos (kz/k)  with  increments  of 
kz/k  equal  to  0.0125  (these  increments  automatically  come 
out  of  the  FFT  when  the  sample  spacing  in  z  is  k/2  and  no 
zero  padding  is  applied).  The  vertical  solid  lines  in  Fig.  9  are 
the  limiting  angles  at  which  the  line  that  begins  at  the  point 
source  and  ends  at  the  far-field  observation  point  touches  the 
upper  and  lower  truncation  edges  of  the  scan  cylinder.  As 
expected,  outside  these  limiting  angles  the  far-field  error  is 
so  large  that  the  computed  far  field  is  useless. 

The  far-field  error  curves  have  spikes  that  correspond  to 
nulls  in  the  probe  patterns  (see  Fig.  3).  These  spikes  are  not 
exactly  symmetric  because  the  point  source  is  not  in  the 
x-y  plane.  They  result  from  the  fact  that  the  probe  output  is 
very  small  at  values  of  z  for  which  the  point  source  is  near 
the  direction  of  a  null.  When  the  far  field  is  computed  in  the 
corresponding  6  direction,  the  correction  factor  is  large  and 
thus  the  errors  caused  by  z  truncation  lead  to  an  inaccurate 
value  for  the  far  field. 

We  now  compute  the  far  field  of  the  point  source  at  0 
=  7t/2  from  measurements  over  the  4>-z  truncated  scanning 
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FIG.  9.  Far-field  error  at  <f>= 0°  computed  from  a  z-truncated  scan  with  <z=30X,  0<  0<36O°,  and  -40X  <z  <39.5X.  The  far-field  error  is  calculated  at  6 
=arccos (kzlk)  with  increments  of  kJk  equal  to  0.0125.  The  limiting  angles  in  9  are  indicated  with  vertical  solid  lines.  The  field  is  produced  by  a  single  point 
source  at  (*,y,z)=(0, 12X,-5X). 


area  given  by  a=30\,  -40\<z<39.5\,  and  60 °<<f> 
<300°.  Figure  10  shows  this  scanning  area  as  a  solid  circu¬ 
lar  arc.  We  shall  compute  the  far  field  with  two  different 
angular  sampling  spacings:  (i)  A<£=1°  resulting  in  A^=360, 


and  (ii)  A  <£=2°  resulting  in  N^=  180.  Hence,  the  relevant 
taper  functions  are  given  by  the  summation  in  Eq.  (37)  trun¬ 
cated  at  |m |  =  180  and  \m\  =90.  Figure  11  shows  upper  enve¬ 
lopes  of  these  taper  functions  for  piston  probes  with  h 


FIG.  10.  Truncated  scan  cylinder  with  radius  a- 30X. 
The  scanned  region  is  the  solid  circular  arc  given  by 
60°  <0<3OO°.  The  point  source  is  located  at  (x,y,z) 
=(0, 12X,-5X),  and  the  two  limiting  angles  are  marked 
by  *. 
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